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We are developing a general, unified, and rigorous analytical framework for using gravitational 
lensing by compact objects to test different theories of gravity beyond the weak-deflection limit. 
In this paper we present the formalism for computing corrections to lensing observables for static, 
spherically symmetric gravity theories in which the corrections to the weak-deflection limit can be 
expanded as a Taylor series in one parameter, namely the gravitational radius of the lens object. 
We take care to derive coordinate-independent expressions and compute quantities that are directly 
observable. We compute series expansions for the observables that are accurate to second-order in 
the ratio e = 'd,/'&E oi the angle subtended by the lens's gravitational radius to the weak-deflection 

1/2 

Einstein radius, which scales with mass as e oc M, . The positions, magnifications, and time delays 
of the individual images have corrections at both first- and second-order in e, as does the differential 
time delay between the two images. Interestingly, we find that the first-order corrections to the total 
magnification and centroid position vanish in all gravity theories that agree with general relativity 
in the weak-deflection limit, but they can remain nonzero in modifled theories that disagree with 
general relativity in the weak-deflection limit. For the Reissner- Nordstrom metric and a related 
metric from heterotic string theory, our formalism reveals an intriguing connection between lensing 
observables and the condition for having a naked singularity, which could provide an observational 
method for testing the existence of such objects. We apply our formalism to the Galactic black 
hole and predict that the corrections to the image positions are at the level of 10 micro- arcseconds, 
while the correction to the time delay is a few hundredths of a second. These corrections would be 
measurable today if a pulsar were found to be lensed by the Galactic black hole; and they should 
be readily detectable with planned missions like MAXIM. 
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I. INTRODUCTION 

The gravitational deflection of light provided one of 
the first observational tests of general relativity. Now it 
is routinely observed in a broad array of astrophysical 
contexts ranging from stars through galaxies and clus- 
ters of galaxies up to the large-scale structure of the uni- 
verse. (See [lyja] for thorough discussions of gravitational 
lensing, and Q for a recent review of astrophysical ap- 
plications.) All of the effects seen so far occur in the 
weak-deflection, quasi-Newtonian limit of general relativ- 
ity. Thus, while gravitational lensing has proven valuable 
for measuring masses in astrophysics, it has not yet been 
able to test fundamental theories of gravity in physics. 

There has been significant theoretical effort over 
several decades to understand lensing in the strong- 
deflection regime (e.g., [Illlllll3,|ili|l3lll|ll), 

in particular for situations in which (i) the lens is com- 
pact, static, and spherically symmetric, (ii) the metric is 
asymptotically flat far from the lens, and (iii) the source 
and observer lie in the asymptotically flat regime. (See 
llOj l for an approach that does not require asymptotic 



flatness.) The resulting theory has yielded a remarkable 
prediction: There should be an infinite series of images 
very close to and on either side of the black hole's pho- 
ton sphere, corresponding to light rays that loop around 
the black hole once, twice, etc. before traveling to the 
observer (e.g., 0, S 0, Sll IHI). The possibility that 
these "relativistic images" could be used to test strong- 
deflection gravity has created considerable excitement, 
inspiring extensive studies of their properties in various 
famihar metrics (e.g., H H i, 0, i, III El) 

as well as 

those arising from string theory and braneworld gravity 
[13L IT3 . Iisj , together with assessments of prospects for 
detecting the images jl^ . 

Unfortunately, the relativistic images are exceedingly 
faint (a flux correction of order 10^^"* for the Galactic 
black hole '3, T^), which makes it important to consider 
whether there are any other observable effects that can 
be used to test theories of gravity. The primary and sec- 
ondary lensed images — which travel from the source 
to observer without looping around the lens — do not 
usually pass close enough to the lens's photon sphere 
to experience extreme strong-deflection lensing, but they 
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may nevertheless be affected by various orders of post- 
post-Newtonian (PPN) correction terms. Since these im- 
ages are much easier to detect, it is valuable to compute 
the corrections to their observable properties. The first- 
order corrections have been studied for the Schwarzschild 
metric (which is standard point-mass lensing the 
Reissner-Nordstrom metric and metrics with gen- 
eral PPN terms 18, 1^ and mass currents |^]. Cer- 
tain aspects of higher-order PPN corrections have been 
studied The lowest-order (i.e., weak-deflection) the- 
ory has been studied for metrics from string theory and 
braneworld gravity [TsL llsl |23 | . Weyl gravity has been 
investigated extensively [2J| and so will not be treated 
here. Unfortunately, much of the existing work seems to 
comprise a diverse collection of results, some of them re- 
discovered several times, that makes it difficult to draw 
general conclusions. Also, most of the previous studies 
have focused on the light bending angle rather than on 
quantities that can be observed directly in extra-solar 
lensing scenarios. 

Our goal is to develop a general framework for com- 
puting corrections to a core set of observable properties 
of the primary and secondary lensed images in a general 
geometric theory of gravity. In the process, we shall unify 
the previous lensing analyses into a common framework, 
and extend them to a higher order of approximation. At 
the same time, we shall demonstrate how to handle gen- 
eral gravity theories using PPN terms up to third-order. 
One crucial part of our formalism is that we take care to 
work with observable quantities and avoid coordinate de- 
pendence. For instance, most previous studies of correc- 
tions to the weak-deflection bending angle have expressed 
the results in terms of the light ray's radial distance of 
closest approach to the lens. However, such a radial dis- 
tance is a coordinate-dependent quantity. Ambiguities 
created by different choices of coordinates can be allevi- 
ated by working instead with the impact parameter of 
the light ray |23|. The impact parameter is an invari- 
ant of the light ray (a constant of motion), and is given 
geometrically by the perpendicular distance, relative to 
initial observers at infinity, from the center of the lens to 
the asymptotic tangent line to the light ray trajectory at 
the observer. It is the quantity to use when defining the 
observable angular position of a lensed image and will 
play a key role in our formalism. 

There are several reasons to compute corrections to a 
higher order than has been done before. One is that cer- 
tain lensing observables have first-order corrections that 
vanish in general relativity (see, e.g., Our higher- 

order formalism will allow us to understand the fortu- 
itous cancellations that cause those terms to vanish, and 
to find the lowest-order non-vanishing corrections. By 
being very general, we will also determine whether the 
cancellations are generic or restricted to specific fami- 
lies of gravity theories; we shall show, for example, that 
the cancellations also occur in all theories of gravity that 
agree with general relativity in the weak-defiection limit. 
A second reason for working to high order is that car- 



rying the expansion far enough will, in principle, allow 
us to bridge the gap that now exists between weak- and 
strong-defiection analyses of lensing by black holes and 
compact objects. 

Our generality will also allow us to find an unexpected 
connection between lensing observables and certain kinds 
of naked singularities. The Reissner-Nordstrom metric is 
usually taken to describe a charged black hole in general 
relativity, but if the charge parameter exceeds a thresh- 
old value, then the metric describes a naked singularity 
instead. We shall discover that the corrections to cer- 
tain lensing observables can have a negative value only 
if the condition for a naked singularity is satisfied. This 
could provide a direct observational test for these exotic 
singularities, which are conventionally ruled out by the 
still-unproved Cosmic Censorship conjecture [25l |. 

This is the first in a series of papers intended to develop 
the complete lensing framework. Here we begin with a 
thorough analysis of lensing by a static, spherically sym- 
metric compact body, using a formalism that can handle 
all gravity theories in which the metric can be expressed 
as a series expansion in the single parameter m, (the 
gravitational radius of the compact body). In separate 
papers we will generalize to metrics with two parame- 
ters and metrics that describe a rotating compact body. 
In addition to presenting the formalism, we will discuss 
some possible astrophysical applications. It is our hope 
that these studies will be relevant to the next generation 
of black hole imagers, such as MAXIM 



II. BASIC ASSUMPTIONS OF THE 
FORMALISM 

Since we seek to lay out a general framework, we should 
begin by stating the assumptions very clearly. Consider 
a compact body of mass M,, perhaps a black hole or 
neutron star, that is described by a geometric theory of 
gravity. This means that the body's gravitational field 
is determined by a spacetime metric appropriate to the 
theory of gravity in question. (See ijlll Al for the form of 
the metric.) Possibilities include general relativity, whose 
spacetime metric obeys the Einstein equation, and mod- 
ified gravity theories, so named because their metrics are 
not governed by Einstein's equation. 

We study how the body acts as a gravitational lens 
by considering light rays that travel directly from the 
source to the observer without looping around the com- 
pact body. Figure ^ gives a schematic diagram of the 
lensing situation and defines standard quantities: B is the 
angular position of the unlensed source; is the angular 
position of an image; d is the bending angle of the light 
ray; and d^, ds, and cIls are the observer-lens, observer- 
source, and lens-source angular diameter distances, re- 
spectively. Also important is the impact parameter 5, 
which is an invariant of the light ray (a constant of mo- 
tion) and is given geometrically by the perpendicular dis- 
tance, relative to inertial observers at infinity, from the 
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FIG. 1: Schematic diagram of the lensing geometry. 



center of the compact body to the asymptotic tangent 
Une to the hght ray trajectory at the observer. As noted 
in ^ when defining observables it is crucial to use the 
invariant impact parameter rather than the coordinate- 
dependent distance of closest approach (also see psj). 
The angular image position is defined in terms of the 
impact parameter as •& = sm~^{b/dL)- 

From the figure, elementary trigonometry establishes 
the relationship (see 

tanB = tan?9 - D (tani? + tan(d - ■&)) . (1) 

where D = dLs/ds- This equation agrees very well with 
the full relativistic formalism for light propagation 9] . In 
this paper, we take it as the general form of the gravita- 
tional lens equation. 

Henceforth, the angles describing image positions are 
assumed to be positive, which then forces the source's 
angular position to take on a positive or negative value 
depending on the image's location. Explicitly, the source 
is assumed to be fixed at an angle B, which is taken to be 
positive when studying an image on the same side of the 
compact object as the source, and negative when study- 
ing an image on the opposite side. (Note the different 
conventions for the signs of i) and B.) 

We make the following assumptions: 

Al: The gravitational lens is compact, static, and 
spherically symmetric, with an asymptotically flat 
spacetime geometry sufhciently far from the lens. 
The spacetime is vacuum outside the lens and flat 
in the absence of the lens. 

A2: The observer and source lie in the asymptotically 
flat regime of the spacetime. 

A3: The light ray's distance of closest approach tq and 
impact parameter b both lie well outside the gravi- 
tational radius m, = GM,/c^, namely, m,/ro <C 1 
and m,/fe ^ 1. The bending angle can then be ex- 
pressed as a series expansion as follows in the single 
quantity m,/6: 




(2) 



The coefficients Ai are independent of m,/6, but 
may include other fixed parameters of the space- 
time. Since b and m, are invariants of the light ray, 
eq. ||2J) is independent of coordinates. Note that 
the subscript of Ai conveniently indicates that the 
component is affiliated with a term of order i in 
m,/6. 

Assumption Al is a natural first step to take for the 
study of lensing by neutron stars and black holes. The 
spherical symmetry conveniently allows us to restrict at- 
tention to light rays moving in a plane. As we shall see 
in ijVIl assumption A2 certainly holds in interesting as- 
trophysical settings. 

As for assumption A3, the form of the general bending 
angle a in eq. Q is a series expansion expressing correc- 
tions to the standard weak-deflection bending angle in 
general relativity, 

-w-(t)-(t)- 

where = tan^^(m,/dL) is the angle subtended by the 
gravitational radius of the compact object, and "w" sim- 
ply indicates that in the last expression we have used 
the standard small-angle approximation. (This approx- 
imation is not valid when we work beyond linear order, 
as discussed in ^IIV Al ) By writing the cofRcient of the 
linear term as Ai , we allow for the possibility that grav- 
ity theories can differ slightly from general relativity in 
the quasi-Newtonian regime. In practice, though, obser- 
vational analyses imply that Ai is quite close to 4: for 
example, recent compilations of 20 years' worth of data 
imply that Ai = 3.99966 ± 0.00090 There are as yet 
no good observational constraints on the higher-order co- 
efficients A2 and A3 in so they are to be determined 
from the gravity theory in question (as discussed below). 

III. LIGHT BENDING ANGLE IN VARIOUS 
GRAVITY THEORIES 

A key part of any lensing framework is the light bend- 
ing angle. We state the general form of the metric we 
consider, and derive the exact, general expression for the 
bending angle. We obtain the desired series expansion 
of the bending angle for the Schwarzschild metric, work- 
ing to higher order than has been done previously, and 
taking care to express the result in terms of the invari- 
ant impact parameter. We then repeat the analysis for 
a general post-post-Newtonian (PPN) metric, and illus- 
trate how various gravity theories can be studied with 
this approach. 

A. Form of the Metric and Bending Angle 

From assumption Al, the spacetime geometry is pos- 
tulated to be static, spatially spherically symmetric, and 
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asymptotically Minkowski: 

ds^ = ~A{f) + B{f) df2 + C{f) dV? , (4) 

where dO^ is the standard unit sphere metric, and 
A{f) 1, B{f) 1, C{f) — > as r ^ oo. Again 
by Al, the metric is Minkowski in the absence of the lens 
and we make the natural mathematical assumptions (for 
regularity in coordinate changes and constancy of causal 
structure) that A, B, C, dA/df, and dC/df are all pos- 
itive in the region outside the lens through which the 
light rays of interest propagate. Note that we do not re- 
quire that the metric approaches the GR weak-deflection 
limit for f sufficiently large. Since dC j df > 0, the func- 
tion C(r) is invertible and allows a new radial coordinate 
This transforms the metric as follows: 

-A{r) di2 + B{r) dr^ + df)^ . (5) 

Because of the spherical symmetry, the geodesies of 
(|5]l lie in a plane, which we can take to be the equatorial 
plane. It then suffices to work with a metric of the form 

ds^ = -A{r)d£- + B{r)Ar^ + d^^ (6) 

where (/? is the azimuthal angle. The Lagrangian / de- 
termining the null geodesies of the metric Q is given 
by 



d.s2 



(7) 

where (wt, w^, Wi^) is a velocity vector with components 
along the 9t, 9^, and coordinate directions, respec- 
tively. Note that this Lagrangian is independent of t 
and ip. The stationary curves are governed by the Euler- 
Lagrange equations, which yield the geodesic equations 
and provide two constants E and L of their motion: 



di ' dip ' 



which give 



Ks) 



E 



0{s) = 



(8) 



(9) 



A{r(s)) ' ' r2(s) ' 

Here, a dot denotes differentiation with respect to an 
affine parameter s along the geodesic. In addition, re- 
stricting the Lagrangian to the position-velocity path of 
the null geodesic (in the tangent bundle over spacetime) 
gives a third constant of the motion: 

- /(r(s),i(s),f(,s),^(s)) (10) 

= -A{r{s))P{s) + B{r{s))r'^{s) + r"^ [s) P^"" {s) . 

This yields 



where h = \L/ E\ and is called the impact parameter. If 
the light ray's distance rp of closest approach to the lens 
occurs at affine parameter s = sq such that r(so) — 0, 
then h and ro are simply related by 

1 _ A(ro) 
62 rl ■ 



(12) 



Equations Q and lfTT|l then yield 



dtp 

IF 



± ■ 



AB 



1/62 _ ^/^2 



(13) 



(since C > and dip/dr — tp/r — ±\(p/r\). The plus 
(minus) sign corresponds to the portion of the light's tra- 
jectory where r increases (decreases) as a function of ip. 
Similarly, equations Q and ((TT|l also give 



dt 
dr 



= ± 



AB 



1 



bA \ - A/r 



(14) 



where the plus (minus) sign corresponds to the portion 
of the orbit where r increases (decreases) as a function 
oft. 

Consider a light ray that originates in the asymptoti- 
cally flat region of spacetime and is deflected by the com- 
pact body before arriving at an observer in the flat region. 
Equation (|13|l then yields the following expression for the 
bending angle (e.g., f28|): 
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2 Y 1/62 - A/r^ 



dr — TT . (15) 



B. General Relativity: Schwarzschild Metric 

1. Schwarzschild Metric 

For a spherical, electrically neutral compact body, gen- 
eral relativity yields a unique metric that is asymptoti- 
cally flat and obeys the vacuum Einstein equation: the 
Schwarzschild metric. It is characterized by the single 
parameter m, and can be written in the form of (jSJ with 



■ , N 2m, , , / 2m, \ , , 

Air) = 1 , B{r) =1 ] . 16 

r \ r 



We need to relate the distance of closest approach tq 
to the impact parameter 6. Using eq. H12() . we can write 
6 in terms of tq: 



6 



1 - 2 



m. 



(17) 



Inverting this to find t-q in terms of 6 yields (cf. [23, p. 
145) 



V3 



33/2 m. 



(18) 
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Since m,/b ^ 1 by A3, we can Taylor expand the right-hand side in powers of m,/6 to obtain 



m, 3 /m. 



m. 



105 /m.^4 



48 



m. 



3003 /m. 



16 



O 



m. 



(19) 



2. Schwarzschild Bending Angle 



From eqs. (|15ll and (|l(j|l . the exact form of the hght bending angle can be written as 

a{ro) = 2 / - TT. (20) 

We change variables to x = tq/w, and also substitute for the impact parameter b using H17|) . This yields 

Q(ro) ^2 f \ '^'^ TT, (21) 

where h = m,/ro. (Note that for computational purposes, a numerically stable expression can be obtained by further 
changing variables to x = cosr/.) For pedagogical purposes, rewrite the integral as 

dx dx ^22) 



Vl-2h- x^ + 2hx^ Jo Vl-x2 ^1 - 2h{l - a;3)/(l - x^) ' 

By assumption A3 we have h < 1/3, which means physically that the light ray is outside the photon sphere at 3m,. 
The rational function (1 — x'^)/{l — x'^) is monotonically increasing on [0, 1] with a maximum value of 3/2 at x = 1. 
This means that < 2h{l — x'^)/{l — x^) < 1, so we can Taylor expand [1 — 2h{l — x'^)/{l — x'^)]~^^^ in a geometric 
series. Carrying out the integration in cq. H22I) term by term then gives 

f . 15 \ , 2 /122 15 \ ^3 / 3465 \ 4 /7783 3465 \ 5 

( 21397 310695 \ ,« .^.,^7 
+ (-^ + ^v^j.« + C.(.)^ (23) 

To obtain the bending angle in terms of the invariant impact parameter 6, we use eq. (|19f) to write tq in terms of h. 
This yields 

where 

1571- , 128 , 34657r , 3584 , 2552557r 
^1=4, A2 = — , ^3 = —, ^4 = ^, ^^-^5^- (2^) 

It follows that eq. (|24() has the form required by assumption A3. We shall not actually use terms beyond O (m,/6)'^, 
but we have included some higher order terms because they have not appeared in the literature before. 



C. PPN Approach 



1. PPN Metric to Third Order 



The post-post-Newtonian (PPN) formalism is a convenient way to handle the wide range of gravity theories in 
which the weak-deflection limit can be expressed as a series expansion in the single variable m,. The formalism was 
extended to third-order in plj . However, we reformulate that treatment to have a better fit with our approach, and 
present new results. 
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(28) 



Express the coefBcients of the standard metric (0) in a PPN series to third-order as follows: 

A(r) = 1 + 2«i + 2a, + 2a, + ... , (26) 

fi(.-)^l-2t,(l) +4!.(1)'-8.,(1)\.... (27) 

where is the three-dimensional Newtonian potential with 

r' 

If the metric is in isotropic form, namely, 

As^ = -A'{r')dt^ + S'(r')[dr'2+r'2df]2], (29) 
then the PPN convention is to write 

AV)^l + 2a'^ + 2f<'(^)%|r(^)\..., (30) 

«V)^.-2y^^.|.(^)'-i,/(^)\..., 

where 4)' /c? — — m,/r'^, and (a', /?', 7', 5', 77') denote the Eddington-Robertson parameters (with primes added to 
avoid confusion with standard lensing quantities). The parameters are chosen so that the Schwarzschild metric has 

a' = /3' = 7' = 5' = e' = r/ = 1 . (32) 

We can relate the parameters in the standard and isotropic forms of the metric by comparing the time, radial, and 
angular parts of the two metrics, yielding the relations 

A{r) = A'{r') , = B'{r') r'^ , B{r) dr^ = B'{r') dr'^ . (33) 

The second and third relations yield 



In r' = / xEIll dr + const , (34) 



where the constant is chosen so that r' /r 1 as ?' 00. Plugging r' into 1)33(1 and identifying terms on the left- and 
right-hand sides that have the same order in m,/r, we find the following correspondence between the standard and 
isotropic coefficients: 

fli — a' , (35) 
bi - 7', (36) 
a2 = /3'-aV, (37) 

62 = ^ , (38) 

3C + 3a'S' ~ 8/3Y -t- 2a^7^^ 
as = ^ , (39) 

+ isyy - 2y3 

63 = ^ • (40) 

The correspondences (I39|l and H4UI) have not been worked out before in the literature. For reference, the Schwarzschild 
metric has 

ai = 61 = 62 = 63 = 1 , a2 = 03 = . (41) 

We need to relate the distance of closest approach ro to the invariant impact parameter b. Using eq. H12() . we can 
write b in terms of tq: 



b = ro 



^ ^ m. ^ 3al - 2q2 / ^ 5aj - 6aiQ2 -I- 203 / rri,\^ Q ( ^ 



(42) 
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To invert this relation and find tq in terms of b, we postulate a relation of the form 



To = b 



m. 



m. 



m. 



m. 



(43) 



plug this into H42|l . and solve for the constants Ci by requiring that the coefficient of each power of m, vanishes. This 
yields 



ro = b 

The result (|44|l is also new. 



1 - ai — 





2a2 — 3af / m. 



(as - 4aia2 + 40^^ 



3^ /m. 



m. 



(44) 



2. PPN 



Annie 



To compute the light bending angle, we take the exact expression H15|) . plug in the PPN metric functions (|26|) - (|27|l . 
and change integration variables to x — ro/w. We also substitute for the impact parameter b using eq. H42|l . Finally, 
we introduce h = m,/ro, so that the series expansion becomes a Taylor series in h. 

The first step yields 



a{h) ^ 2 



dx 



l + h{bix 



1 



1 + x 



3al 



^ 2^ ^ ^' 2(l + x)2 



+ {2al + aibi) 



X 



1 + x 



2 



ai{8al + Aaibi ~ bf + Ab2)ix - 1) + 4aia2(l - 2x) + 2x{a3 - 0261) 



-(&? -46162 + 863)^^ + 



5af 3af(4ai+6i) 2a3 



(l + a;)3 (1 + ^)2 



1 + x 



ai(20af - 10a2 + 7ai6i - 6f + 462 
1 + a; 



(45) 



Carrying out the integration term by term then yields 



a{h) = 2(ai +6i)/i + 
"67 



2ai(ai + 61) + ( 2al — 02 + ai6i 7 + 62 ) tt 



^af - 18aia2 + 4a3 + 9a?6i - 26i(a2 + ai6i) + 80162 + 3 " 46i62 



-ai( 4ai - 2a2 + 2ai6i - y + 262 )7r 



(46) 



which agrees with [23|. However, the expression H46|l is coordinate-dependent since it is written in terms of the 
distance of closest approach. To obtain an invariant expression (so we can discuss observable quantities), we use (|44f) 
to write ro in terms of 6 and obtain: 



d(6) = A, 



where 



A, 



Ai = 2(ai-)-6i), 

A2 = ( 2a? - 02 +ai6i - ^ -^62 ) TT, 



-It 



35a? + 15a?6i - 3ai (l0a2 + bj ~ 462) -I- 603 + bf - 60261 - 46162 



(47) 

(48) 
(49) 

(50) 



This invariant bending angle expression has not appeared in the literature before. 
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D. Sample Gravity Theories 



In general relativity, the Reissner-Nordstrom metric describes a black hole with physical charge Q. Although charge 
is generally expected to become neutralized, in certain circumstances it could persist for 1,000 to 10,000 years |30| . 
perhaps allowing detection of such black holes. The Reissner-Nordstrom metric has the form of ([SJ with the metric 
functions 

A{r)=B{r)-' = 1-—^ + ^ , (51) 

where q = VG Q/{cPm,) is a dimensionless parameter. If q^ > 1, the metric actually describes a naked singularity 
rather than a charged black hole (e.g., [l3ill3)- It is natural to view H51|) as a series expansion in m,/r. We can then 
identify the PPN coefficients: 

9 9 2 

q q q 

ai = l, 02 = y, 03 = 0, 6i = l, 62 = 1- ^, b3 = l- — . (52) 

Using eqs. ()48 |I -H50 () . we can determine the coefficients in the expansion of the bending angle: 

Ai=4, A, = i5-q')^, A,^^-16q\ (53) 

Notice that if 5 ^ we recover the results for the Schwarschild metric. Also notice that the condition for a naked 
singularity corresponds to A2 < Stt. Turning this around, we may say that for the Reissner-Nordstrom metric of 
general relativity, the coefficient A2 can he negative only if there is a naked singularity. 

In heterotic string theory, modifications of the Einstein equation lead to a different charged black hole solution, 
which is often called the Gibbons-Maeda-Garfinkle-Horowitz-Strominger (GMGHS) black hole ^SlJ. In the 
metric's lensing properties were computed to lowest order in both the weak- and strong-deflection regimes, although 
only with a coordinate-dependent approach. We can include the weak-deflection limit in our invariant formalism, and 
thereby obtain new results. For a black hole with gravitational radius m, and charge Q, the GMGHS metric is often 
written in the form 



d.^ = - (^1 - ^) dt- + (1 - ^) " + (1 - ^ j dO^ , (54) 

where again q = ^/GQ / {c?m,). As in the Reissner-Nordstrom case, this metric describes a naked singularity if the 
charge parameter exceeds some threshold, in this case q^ > 2. We convert to standard coordinates by setting 

2 _-2 f. 



r'^f'il- ^ . (55) 



This puts the metric into the form of (jSJ with the metric functions 



A{r) = 1 - ^ , ^ , (56) 



B{r) = 1 - 



q^m, + \J q'^mi + 4r^ 

4m, \ 4r^ 



q^m, + yjq^ml + 4r'^ I mlq"^ + 4r'^ 



(57) 



Expanding the metric functions as Taylor series in m»/r, we can identify the PPN parameters 



2 4 2 4 2 4 

ai = l, 02 = -, 03 = -, h = l, 62 = 1----, b,^l^---. (58) 

The coefficients in the expansion of the bending angle are then 

4i=4, A2 = {60- Uq'-q^)^, A, = ^-16q\ (59) 

In this case, the condition for a naked singularity corresponds to A2 < 2tt. Or, we may again say that for the GMGHS 
metric A2 can be negative only if there is a naked singularity. 



9 



IV. LENSING FRAMEWORK IN VARIOUS GRAVITY THEORIES 

We can now move beyond the bending angle to compute corrections to the observable properties of the primary 
and secondary lensed images. In this section we focus on the positions and magnifications of the images; time delays 
are deferred to ^ 

A. Lens Equation 

We start with the general lens equation, 

tanS = tan?9 - D (tani? + tan(d - z?)), (60) 

and seek an appropriate series expansion. First, we change variables to match the scalings commonly used in the 
astrophysical lensing literature. A natural scale is the weak-deflection angular Einstein ring radius. 



AGM.dLS 



We then define: 



4D 



(61) 



(62) 



In other words, /3 and 9 are the scaled angular positions of the source and image, respectively. The quantity e 
represents the angle subtended by the gravitational radius normalized by the angular Einstein radius, and it becomes 
our new expansion parameter. 

The second step is to postulate that the solution of the lens equation can be written as a series expansion of the 
form 



= oa 



(63) 



where is expected to be the image position in the weak-deflection limit, and the coefficients 9i and 02 of the 
correction terms remain to be determined. (This is a standard perturbation theory analysis; e.g., (32i] . ] 
After making these substitutions, we flrst find that we can write the bending angle as 



Ai A2-Aiei 2 



1 

"^0 



A3-2A291+A1 -D 



(64) 



Note that since we are expanding beyond linear order, it is important to use the exact geometric relations between 
physical and angular radii: i9 — sin~^ (b/dL) and i?, = tan^^(m,/(iL). The standard small angle approximations 
(i? Ri b/dL and « m./d/,) are valid only at linear order. 
Now making the substitutions in the lens equation and Taylor expanding in e, we find: 







D 



-if3 + 400 



D 



-A2 + (Ai + Ael) Oi 



-^\aI + 3A3 - 12^2^612 + Ai(56L>26i4 + Wl - 30o02) - 2 32D'^e'f,{0f, - /?^) + 3A26I1 + 69p2 



(65) 



This is the desired series expansion of the lens equation. 



B. Image Positions 



We now solve eq. (|65|) term by term to find the coefficients 9i in the series expansion for the image position. The 
idea is to fix the source position /3 and flnd the values of 9i that make each term in ()65fl vanish. The first-order term 
is just the standard weak-defiection lens equation. 







-(3- 



1 

To ' 



(66) 
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which yields the weak-deflection image position 



(67) 



We neglect the negative solution because we have explicitly specified that angles describing image positions are 
positive. We find the positive-parity image 9q , which lies on the same side of the lens as the source, by using a 
positive angular source position, /? > 0. We then find the negative-parity image 9q , which lies on the opposite side of 
the lens from the source, by using /3 < 0. In other words, we can rewrite H67|l as 



^0 



4 + /32 ± 



One curious feature of the positive- and negative-parity image positions in the weak-deflection limit is that 



^0 '^o 



1. 



Note that 1)68(1 makes it clear that H69|) does not depend on the sign of /3. 

Next, we choose 9i to make the second-order term in ((65(1 vanish. This yields 



A. 



(68) 



(69) 



(70) 



For the third-order term, we still have given by H67|l and 9i given by ((7()|l . and we must choose 62 to make the 
coefficient of in (|65() vanish. This yields 



1 



Ai{A\ - iAl + - 64AiL>2) 

&Al - 2Ai{Al + 3A3) + iA\D - 16AiL>2(3Ai - 8) 
2Al{l - 6D) + 6A3 - 1281)2 + AiD^{192 + Ai(-24 -I- 7Ai)) 



3 00 (Ai +4(?2)3 
- 4 



64D 



48D - 24AiD + A?(-3 -I- 7D) 



128D- 



'(-24 + 7Ai)08|. 



(71) 



Note that we have used the relation f3 — Oq — I/Oq to replace /3 on the right-hand side. It is possible to rewrite the 
right-hand side in terms of (3 instead of 9o, but at this point we believe that is less useful. 

For gravity theories with Ai = A (including general relativity) the correction terms simplify slightly to 



A. 



Ail + 81) ' 
1 



3A^(1 + 29'^) + 4(1 + 0af[ 3A3 + 64 - 64D'^ - 192D0'^ + 192D'^0'^ + i2D^e\ 



It is worth pointing out that the image position can be written to first-order in e as 



-,2/)2 



-i2/i4 



Ao 



Ai 



Aei 



0{eY 



(72) 
(73) 

(74) 



The sign of the first-order correction to the weak-deflection image position is given by the sign of A2 (for both the 
positive- and negative-parity images) . We have already seen that the sign of A2 is connected with the possible presence 
of naked singularities: in the Reissner-Nordstrom and GMGHS metrics, A2 can be negative only if there is a naked 
singularity (see i illl Dp . Thus, the correction to the image positions can provide a possible observational test for naked 
singularities. 



C. Magnifications 



The signed magnification /i of a lensed image at angular position is 



sinS(i?) dB{§) 



(75) 
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We change to our scaled variables as in eqs. 1)62(1 and 1)63(1 . make a Taylor expansion in e, and then substitute for 6i 
and 02 using l)7U() and 1)71(1 . In this way we obtain a series expansion for the magnification: 



where 

fJ-i = 
fJ-2 = 



= Mo + Ml e + M2 + O (e) , 



ml 



16^2613 



A\D^ {bl2-QA{) + AAief^ -384AiD^+Af(4+12£i-9£i^)+ 4(3^3 



3(^1-402)2(^^+4^,2)5 ^ 

+ 64^16';^ \a{{1 + SD) - 9Al + 3^11^3 + li^AiD"^ - 192AlD^ + 24:AID^ 
'9Al + 3A1A3 + 96AId^ - 48AId^ + A^{1 + 3D + 2D^) 



(76) 

(77) 
(78) 
256i:'2 



- 256 6*^ 



256 6ig 



12^3 - 96AiD^ + Ai{A + 12D + 7D^) 



l02AAlD^el^ 



(79) 



Recall that the sign of the signed magnification indicates the parity of the image: /j, > for 9^ , the positive-parity 
primary image; while /i < for 9^ , the negative-parity secondary image. 
For a gravity theory with Ai — 4, these expressions reduce to 



Mo 



Ml 



M2 



0^ 
''o 



A29'o 



4(1 + 02) 



2\3 



24(02 _ i)(^2 + i)s 



768D9l{9l + 1)2 - 64:D^{9l + if (l + 166'g + 9^ 



+9l 256 -I- (512 - 9Al 



2\n2 



2566^ 



12^3(02 



(80) 
(81) 

(82) 



It is worth pointing out that if Ai = 4 the first-order changes in the magnifications of the positive- and negative- 
parity images have the following relation: 



M^ 



A. 



4[1 



9+)2]3 



4[(0„-)2 + l]3 



Ml 



(83) 



In the second equality we used 1)69(1 to write 9q — 1/9q. In gravity theories with Ai = i and A2 > 0, the ^1 
perturbation is negative, so it makes ^'^ less positive (fainter) and /i^ more negative (brighter) by exactly the same 
amount. (The opposite occurs for Ai = 4 and A2 < 0.) In other words, the magnifications of the positive- and 
negative-parity images are shifted by the same amount but in the opposite sense. This has important implications 
for the total magnification and centroid (see below). 

Notice that the sign of the first-order magnification correction depends on the sign of A2 , so this provides another 
possible observational test for naked singularities. 



D. Total Magnification and Centroid 

If the two images are too close together to be resolved (as in microlensing) , the main observables are the total mag- 
nification and the magnification-weighted centroid position. Using our results above, we find the total magnification 
Mtot = Im^I + Im~I to be 
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The second-order term can be worked out from (|79|l if desired, but is too complicated to write here. Notice that the 
sign of A2 again determines the sign of the first-order correction. 

It is important to see the factor of {Ai — 4) multiplying the first-order term in H84() . It means that the first-order 
correction to the total magnification vanishes with full generality in any gravity theory with Ai = 4. In such theories, 
the first-order changes in the magnifications of the positive- and negative-parity images exactly cancel (see ea. l83f) . A 
correction then enters at second-order, which is given by 



2 + /3 

A*tot 



^ ^2 )5/2 1 9-^2 - 1 2^3 (4 + ) - 64(4 + ) [i + 12D - {18 + (3^)D^]} + O {ef . (85) 



Here we have found it convenient to express the result in terms of the source angular position (3. 
The magnification-weighted centroid position is defined by 

e+\^l+\-9-\^^\ _ 9+^1+ + g-^- 

fccnt — 1 , I , I I — ■ [ob) 

(The sign in the numerator may be understood by recalling that we use positive angles, 9± > 0.) Our perturbation 
analysis then yields 

o _ Aj{0l- 01 + 01 -01 + 1) -160^0 

- Al0o{0',-0l + 0l-l) ^''^ 



A2{Ai~A)9l 



Ai{16 + AAi + Ai){l + 91'') - 4Ai{32 - Al)9l{l + 9^) 



At{9f, 1)(04 + i)2(Ai + 40l){4 + A,0l) 

+ (256 - 64Ai - 32Al + SAl + Aj)0^{l + 0^) + (256 + 32 Al + 2AAl + Al)0^ 



0{ef 



As before the second-order term can be worked out but is too complicated to write here. Notice again that the sign 
of A2 yet again determines the sign of the first-order correction. 

We again see a factor of {Ai — 4) multiplying the first-order term. Thus, the first- order correction to the centroid 
vanishes exactly in any gravity theory with Ai — 4. In such theories, a correction appears at second-order and is given 
in terms of /3 by 



2 + /32 24(4 + /32)(2 -K/32)2 



eeent = ^^^^^ - .... , .2Vo , - 12.43 (4 + ) ^ 128(4 + ) (2 - ) 



-64(4 + 13') [(9 + I3')D- 6] D(3' j e' + O (e)^ . (88) 
Notice that in this case {Ai ~ 4), the sign of A2 does not affect the corrections, but the sign of A3 does. 

V. TIME DELAYS IN VARIOUS GRAVITY THEORIES 

Having studied the positions and brightnesses of the primary and secondary images, we are now ready to compute 
their time delays. The analysis parallels some of what has gone before, but is different enough to warrant a separate 
treatment. 



A. General Formalism 

Let Rsrc and i?obs be the radial coordinates of the source and observer, respectively. From geometry relative to the 
flat metric of the distant observer, who is assumed to be at rest in the natural coordinates of the metric lO, we can 
work out (see Figure ^ 

i?obs = dL , i?src = {dls + dl tan^ B) . (89) 
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The radial distances are very nearly the same as angular diameter distances since the source and observer are in 
the asymptotically flat region of the spacetime. In other words, the distortions in distances near the black hole are 
assumed to have little impact on the total flat metric distance from the compact body to the observer or source. We 
focus on spacetimes that would be flat in the absence of the lens (see assumption Al), and in that case the light ray 
would travel along a linear path from the source to the observer with length ds/ cosB. 

The time delay is the difference between the light travel time for the actual ray, and the travel time for the ray the 
light would have taken had the lens been absent. This can be written as 



CT = T{Rs 



T{Robs) - 



ds 
cosB 



with 



T(i?) = f 



dr = 



A{w)B{w) 



A{w) y l/fe2 - A{w)/i 



■ dw . 



(90) 



(91) 



where we have used eq. (|14f) for dt/dr. Unlike eq. (|15|) . this integral cannot extend to infinity because the travel time 
would diverge. 

Note that there are no redshift factors in these equations because we are assuming a spacetime that is static and 
asymptotically flat. This assumption is not overly restrictive because expected applications involve non-cosmological 
lens systems (see 



B. General Relativity: Schwarzschild Metric 

To compute T{R) for the Schwarzschild metric, we take the metric functions from eq. H16|l . change integration 
variables to x — rt^/w, and use H17|l to relate 6 to rg, so that we find 



T(R) = ro [ (1 _ 2/ix)"^ (1 



l-2hx 
1 - 2/i 



^1/2 



dx , 



where h = m./rp. Expanding the integrand as a Taylor series in h and integrating term by term yields 
T(i?) = Jr^ - tI + hro 



+ h^ro 



(92) 



(93) 



2(1+0' 



15/7r_ . _i A (60 + 157g + 133$^ + 35^^) ^1 - ^ 
2 V2 V + 2(1+^3 



where ^ — tq/R- 

This expression is currently written in terms of the distance of closest approach ro, but we want to rewrite it in 
terms of the invariant impact parameter b. We can use (|19() to make the translation, but we must then consider the 
nature of the series expansion. The expansion in ft, = m,/ro naturally becomes an expansion in m,/5. It may be less 
obvious, but we also want to expand in b/R, which is small because both the source and observer lie far from the 
lens. To do the joint expansion, we need to consider the amplitudes of m,/5 and b/R. In terms of angular variables, 
we have two fairly simple cases. 



^obs 



m, 

T 

b 

dL 







De. 



and one that is slightly more involved, 
b b 



1 -D 



D{l-D) 



Rs 



dls + dlUn^B VD^+tan'B ds ^D^+tan'B dL 



tan^ B 



(94) 
(95) 

(96) 
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The bottom line is that m,/5, b/Rohs, and 6/i?src are all similar in amplitude, up to factors that depend on D. We 
will do the rigorous series expansion of the time delay momentarily, but for now it is instructive to expand (|93|l in 
both m,/6 and b/R, taking care to collect terms of a given order in any combination of the two quantities. Working 
to third-order, we find: 



T{R) 
R 



1 b 

2 R 



b 
R 



m. 



-4^1nfiA 
b \2 R 



+ 



Wtt b 
~ R 



m. 



(97) 



There is no term that is linear in m,/& or b/R. The first term in braces is of second-order, while the second is of 
third-order. We have neglected terms of fourth-order and higher. 

That was meant to be pedagogical. To properly compute the series expansion for the full time delay H90|) . we first 
compute T(i?src) and T(i?obs) by plugging (|5^ into We replace ro with b using ((T^ . We change to angular 

variables using b = di,sinz5, and then reintroduce the scaled angular variables 9 and /3 defined in H62I) . Finally, we 
take a formal Taylor series in our expansion parameter e. Putting the pieces together, we find: 



dL dLS 



Ad 



LS 



157r-8(l + 6i2)6li 



(98) 



The interpretation of this result becomes more clear when we recognize that a characteristic scale for the time delay 
is 



_ dbds ,2 . rn. 

cdLS c 



(99) 



Although it is not obvious from the definition, the second equality shows that te is independent of the distances. 
Using this definition, we can write (|98|l as 



r 

TE 



l + P'^-el- In 



dL Bl-dl 

A dLS 



157r-8(l 



Now rewrite 9i (the first-order correction to the image position) using H70() . with Ai 
Schwarzschild metric, to obtain 



T 
TE 



I + 13'' -el- In 



A dLS 



157r 
16^ 



0{e) 



Hs)\ (100) 
4 and A2 = 157r/4 for the 

(101) 



Though not yet apparent, the zeroth-order term in (|101|l reduces to the familiar lensing time delay in the weak- 
deflection limit of general relativity. Rearranging, we can write the zeroth-order term as 



T = 



dL ds 
cdLS 



(102) 



where we have used the identity f3 — 9o — 9q^ from and we have also used the definitions 6*0 — '&o/'&e and 

/3 — B/i^E from H()2|) . The "contant" term C in this expression is independent of i?o and B; it depends on -dE and the 
distances as C = T£;[1 -t- ln(di/4(iis)]/2. 

Finally, to make contact with conventional calculations we examine the differential time delay At = r_ — t+ between 
the positive-parity primary image and the negative-parity secondary image. This can be written as 



At = Aro + eAn + O {e)\ 



where 



Aro = TE 
An = TE 



(0^ 



16 9+9^ 



-'"'if 



(103) 

(104) 
(105) 



It is possible to derive the second-order correction terms in eqs. (|101|l and (|103|) . but they are more complicated. 
Since we have already found the first correction terms for both absolute and differential time delays, it seems unnec- 
essary to write down the higher-order terms. 



15 



C. Time Delay via the PPN Metric 



To compute T{R) for the PPN approach, we take the metric functions from eqs. I|26|l and (|27|l . use H42|l to relate b 
to rg, change integration variables to a; = tq/w, and finally carry out the integration to find 



T(i?) 



«i + (ai +^i)ln I ^ 



(106) 



4al - 2a2 + 2aibi - + 262 ) ( 77 - sin"^ ^ 



-ai ( ai + &i + -fli^ + 61C 



+ h^ro 
+ 0{h)\ 



2(1 + 0' 



2(1 + 0' 



ai 4a^ - 2a2 + 2ai6i - -bf + 2b: 



sin-^C 



where ^ = tq/R. Here -F(0 is a cubic polynomial in ^, which is long and not particularly enlightening so we have not 
written it out (but we do write complete expressions in what follows). 

As in the Schwarzschild case, it is instructive to convert from the distance of closest approach rg to the impact 
parameter b (now using ea . I44|l . and then expand in both m,/& and b/R to find 



T{R) 
R 



1 b 

2 R 



b „ rn, „ , , , m, , /l 6 

^-2a,- + 2(a, + M-ln -- 



:8a? - 4a2 + 4ai6i - bj + 462) J ;^ ' 



However, again we need to be more careful to obtain a rigorous series expansion of the time delay. Repeating the 
analysis discussed around eqs. H98(l - p()()|l . we obtain 



T _ 1 
^ " 2 



, fl2 /i2 «1 + 01 1 

ai+P — In 





1 - 


V 4dL5 J\ 





{8al~Aa2+4.aibi-bl+Ab2)7:-4.{ai + bi + 2el)9i e + O (e)' . 

(108) 

We can rewrite 9i (the first-order correction to the image position) using (|7()|l . with Ai and A2 given by (|48|l and 
(H^ . to obtain 



r _ 1 
^ " 2 



ai + - eg - n-L—Il In 



ai+bi fdL 9^ t?| 



f 8a? - 4a2 + 4ai6i - bj + 462) e + O {ef . (109) 
lb t'n V >^ 



Finally, the differential time delay between the positive- and negative-parity images is 



At = Aro + eAn + O (e) , 



where 



Aro = TE 



An = TE ^ (Sal - 4a2 + 4ai6i - bf + 462) ^^'j,, /° 
Id V 9q 9q 



(110) 

(111) 
(112) 



VI. APPLICATION TO THE GALACTIC 
BLACK HOLE 

To illustrate our results, we consider gravitational 
lensing by the supermassive black hole at the center 
of our Galaxy. The black hole has a mass of M, = 



(3.6 ± 0.2) X 10^ M0 Isjl and is located at a distance 
of (11 = 7.9 ± 0.4 kpc [33 from Earth. (For illustration 
purposes, we adopt the nominal values and neglect the 
small uncertainties.) Its gravitational radius is therefore 
m, = 1.1 X 10^^ cm — 3.4 x 10^^ pc, which corresponds 
to an angle of 6, = 9.0 /las (micro-arcseconds). 
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Suppose the source lies a distance d^s beyond the 
black hole. Typical distances will be d^s ^ 1-100 pc, so 
we have w ^ c^ls ^ rn» which means that both 
the source and the observer lie in the asymptotically flat 
regime of the spacetime, confirming assumption A2. The 
angular Einstein radius is — 0.068 (dis/lO pc)^/^ as, 
which is much larger than 6,. Since the primary and 
secondary images both lie near the Einstein radius, each 
has a light ray with an impact parameter and distance 
of closest approach that are much larger than the grav- 
itational radius, confirming that part of assumption A3. 
Put another way, our dimcnsionless expansion parameter 
is 

^ = f = 1.3x10- xf-^)"'', (113) 

which is small enough to justify working with series ex- 
pansions in e. The natural lensing time scale is te — 71 
sec (independent of the distances; see eg. I99f) . 

We obviously want to compute the corrections to lens- 
ing observables for the Schwarzschild metric in general 
relativity. For comparison, we also consider three other 
cases: a charged black hole described by the Reissner- 
Nordstrom metric in general relativity, with charge pa- 
rameter q = 0.5; and a naked singularity described by ei- 
ther the Reissner-Nordstrom metric of general relativity 
or the GMGHS metric of heterotic string theory (mod- 
ified gravity), with "charge" parameter q = 2.5. (A 
charged black hole described by the GMGHS metric with 
q = 0.5 is very similar to the Reissner-Nordstrom case 
with q — 0.5, so we do not show it.) Note that all four 
cases are from two gravity theories in which the coeffi- 
cient of the leading term in the expansion of the bending 
angle is Ai = 4. 

Figures |21 and 13 show the corrections to the individ- 
ual image positions and magnifications, as a function of 
the scaled angular source position (3, for the four sample 
cases. The scaled first-order correction term for the im- 
age position 9i is of order unity, so the full correction to 
the image position is of order eSe ^ 10 /ias. The cor- 
rection is small, but detectable with high-resolution ra- 
dio interferometry; and it should be readily measurable 
with planned microarcsecond-resolution missions such as 
MAXIM |2g. While the difference between neutral 
(Schwarzschild) and charged (Reissner-Nordstrom) black 
hole cases is fairly small, the naked singularity case 
stands out for having the opposite sign (see ^IVB|) . In 
other words, it will be challenging but not impossible to 
detect the corrections to the usual weak-deflection lens- 
ing. If the corrections can be detected, high precision 
will be needed to distinguish between different kinds of 
black holes. However, it will be easy to rule out certain 
kinds of naked singularities observationally. 

Figures 0] and |31 show the corrections to the total mag- 
nification and magnification-weighted centroid position. 
The corrections to these observables are much smaller — 
the centroid correction is of order e'^'Oe ~ 10~® as — 
because the first-order terms vanish in gravity theories 




-10 r 

-20 y 



-2-10 1 2 

FIG. 2; Terms in the series expansion 1631 for the angu- 
lar image position, as a function of the angular source po- 
sition. Recall that /3 > corresponds to the positive-parity 
image while /3 < corresponds to the negative-parity im- 
age. Sources very close to the origin < 0.1 are not 
shown. (Top) The zeroth-order image position ^o, which is 
the same for all gravity theories. (Middle) The first-order 
term 9i; recall that the full correction term is e9i. Solid 
curve: Schwarzschild metric in general relativity. Dotted and 
short-dashed curves: Reissner-Nordstrom metric in general 
relativity, with charge parameter q = 0.5 and 2.5, respectively. 
Long-dashed curve: GMGHS metric in string theory (modi- 
fied gravity), with charge parameter q = 2.5. (The GMGHS 
case with q — 0.5 is very similar to the Reissner-Nordstrom 
case with q = 0.5.) (Bottom) The second-order term 82. 
Again, the full correction term is e^92- All angular lengths are 
in units of the angular Einstein radius For the Galactic 
black hole, i9_b = 0.068 (dLs/10 pc)'^''^ as and the dimension- 
less expansion parameter is e = 1.3 x 10~* x (dLs/W pc)^^''^. 



with Ai = 4. Although interesting, these quantities will 
be very challenging to use in the near future for realistic 
observational tests of gravity theories. 

Finally, Figure shows the corrections to the differen- 
tial time delay between the positive- and negative-parity 
images. Remarkably, the first-order correction e Ati is 
predicted to be as large as a few hundredths of a second 
— which would be relatively easy to measure if we could 
find a pulsar lensed by the Galactic black hole. As with 
the image positions, distinguishing between charged and 
neutral black holes would require higher precision (not 
unfeasible in the case of a pulsar source), but ruling out 
certain kinds of naked singularities would be easy be- 
cause they lead to corrections with the opposite sign. 
Although we have not explicitly computed the second- 
order corrections, we can estimate that they would be of 
order e'^te ^ 10~® s. 
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0.5 1 1.5 2 

FIG. 5: Similar to Figure |1J but showing the zeroth- 
and second-order terms in the series expansion (1881 of the 
magnification- weighted centroid position, in units of i9_b. 
Again, for our sample gravity theories having Ai = 4, the 
first-order correction term vanishes identically and is not 
shown. 



FIG. 3; Similar to Figure |21 but showing the terms in the 
series expansion 1761 of the individual image magnifications. 
Again recall that the full first- and second-order correction 
terms are e fii and e^fi2, respectively. 
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FIG. 4: Zeroth- and second-order terms in the series expan- 
sion 18511 of the total magnification. Our examples involve 
gravity theories with Ai — 4, so the first-order correction 
term vanishes identically and is not shown. For each /3 > 0, 
we find the positive-parity image using the source position /3, 
and the negative-parity image using the source position — /3, 
and then combine them to obtain the total magnification. The 
line types are the same as in Figure |5| 



VII. CONCLUSIONS 

We have introduced a rigorous and general framework 
for computing corrections to standard weak-deflection 
lensing observables. In this paper we have presented a 
formalism for handling any static, spherically symmet- 
ric theory of gravity in which corrections to the weak- 
deflection regime can be expanded as a Taylor series in 
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FIG. 6: Zeroth- and first-order terms in the series expansion 
lllUI I of the differential time delay between the two images. 
Note that we have defined Arc and An to have dimensions 
of time, so they are expressed in seconds. Recall that the full 
first-order correction is e Ari , where the expansion parameter 
is e = 1.3x 10"'' X (dis/lO pc)"^'^^ for the Galactic black hole. 

the gravitational radius of the compact body acting as a 
gravitational lens. Conceptually, three points distinguish 
our framework from previous studies. First, we take care 
to avoid coordinate dependence by expressing our results 
in terms of invariant quantities. Second, we go beyond 
the bending angle to study quantities that are directly 
observable for extra-solar lensing studies, including the 
positions, magnifications, and time delays of the lensed 
images. Third, our general approach allows us to unify 
the diverse results that have been presented before. 

Besides the framework itself, our main results are series 
expansions for the lensed image positions and magnifica- 
tions that are accurate to second-order in e = t?, 
or the ratio of the angle subtended by the lens's gravita- 
tional radius to the (weak-deflection) angular Einstein ra- 
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dius, as well as series expansions for the lensing time de- 
lays with corrections at first-order in e. The signs of the 
first-order corrections are determined by the sign of the 
coefficient A2 in the series expansion of the light bending 
angle. If A2 > 0, both the positive- and negative-parity 
images are shifted away from the lens, and the positive- 
parity image gets fainter while the negative-parity image 
gets brighter. If A2 < 0, the opposite occurs. 

The sign of A2 appears to be connected with the possi- 
ble existence of naked singularities. In the sample gravity 
theories we have considered — namely the Schwarzshild 
and Reissner-Nordstrom metrics in general relativity, and 
the GMGHS metric in heterotic string theory — A2 can 
be negative only if naked singularities exist. This con- 
nection between hypothetical naked singularities and ob- 
servable lensing quantities is exciting because it offers 
the possibility that (certain kinds of) naked singularities 
can be ruled out observationally rather than just by the 
still-unproved Cosmic Censorship conjecture. Note that 
we must use the qualifying phrase "certain kinds of" be- 
cause not all naked singularities lead to A2 < 0. In these 
less extreme cases, we must go beyond a gross feature 
like the sign of A2 and consider the extent to which fine, 
quantitative constraints on A2 from realistic data could 
rule out all or at least most interesting types of naked 
singularities. 

It has been known that for the Schwarzschild metric 
the first-order corrections to the total magnification and 
magnification- weighted centroid vanish (see We 
have shown that this result is not generic: it depends 
on precise cancellations that can occur only for gravity 
theories in which the leading-order term in the series ex- 
pansion of the bending angle has the value Ai = 4. In 
practice, Ai does equal 4 in the sample gravity theories 



we have considered, and it is constrained by observational 
data to be quite close to 4 (e.g.,|23|). Nevertheless, it is 
important to understand why these particular first-order 
corrections vanish, and to recognize that it is possible to 
devise gravity theories for which that is not the case. 

We have applied our formalism to lensing by the Galac- 
tic black hole. We predict the corrections to the image 
positions to be at the level of 10 microarcseconds, and 
the correction to the time delay between the images to 
be a few hundredths of a second. The position corrections 
would be measurable today with radio interferometry, if 
we could find a radio source that is lensed by the Galactic 
black hole. The time delay correction would be measur- 
able if the source were a pulsar. Even if such a convenient 
source cannot be found, the corrections should be mea- 
surable with planned missions such as MAXIM [26l| . 
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